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ELASTOHYDRODYNAMIC ANALYSIS USING A POWER LAW 
PRESSURE-VISCOSITY RELATION 
by Stuart H. Loewenthal and Erwin V. Zaretsky 

Lewis Research Center and 
U. S. Army Air Mobility R&D Laboratory 

SUMMARY 

An isothermal elastohydrodynamic (EHD) inlet analysis of the Grubin type which 
considers a power law pressure-viscosity relation (p = p Q (l + Kp) n ) and a finite pres- 
sure at the inlet of the contact zone was performed. Numerical results from this analy- 
sis were compared to EHD film thickness data from X-ray measurements made with a 

synthetic paraffinic oil and to results from conventional EHD analysis. 

g 

In the high-contact stress regime; that is, above 1.04x10 newtons per square 
meter (150 000 psi), minimum film thickness as predicted by both present and conven- 
tional theories was consistently less sensitive to applied load (contact stress) than the 
measured data. The present EHD theory exhibits a slightly stronger load dependence 
than do previous isothermal EHD theories which consider a straight exponential 
pressure-viscosity relation but far less than that exhibited by the measured data. 


INTRODUCTION 

The effectiveness of an elastohydrodynamic (EHD) film lies in its ability to main- 
tain separation between two loaded surfaces in moving contact. The size of this spacing 
or film thickness, as it is more commonly referred to, is primarily dependent on the 
rheological properties of the lubricant, the elastic properties of the contacting surfaces, 
and the conditions of operation. Because of its importance on the design life of contact- 
ing machine elements, prediction of EHD film thickness has been the focal point of 
many theoretical and experimental investigations. A compendium of the advancements 
made in EHD lubrication up to about 1964 appears in a text by Dowson and Higginson 
(ref. 1) and more recently in a report by McGrew et al. (ref. 2). 

The bulk of the experimental work conducted in EHD lubrication has been confined 



to conditions of moderate speeds, that is, up to 25.4 meters per second (1000 in. /sec), 

Q 

and moderate loads, that is, maximum Hertz stresses to 1.24x10 newtons per square 
meter (180 000 psi) (refs. 3 to 6). Recently Parker and Kannel (ref. 7) have extended 
EHD film thickness measurements with an X-ray disk apparatus to a range of speed and 
load conditions, that is, surface speeds to 37.6 meters per second (1480 in. /sec) and 
maximum Hertz stress to 2.42x10^ newtons per square meter (350 000 psi) which are 
more in line with current bearing and gear design requirements. Whereas, the film 
thickness measurements determined at conditions of moderate loads and speeds by pre- 
vious investigators have shown reasonably good agreement with EHD theory, the recent 
X-ray data is at odds with theory. In particular, at high contact pressures the sensi- 
tivity of film thickness to load as determined experimentally is far greater than pre- 
dicted by the theory of either Grubin (ref. 8) or Dowson and Higginson (ref. 9). Conse- 
quently, predicted film thickness as at high loads are substantially higher than those 
observed . 

In an attempt to resolve this discrepancy and to determine whether the X-ray meas- 
surements were themselves in error, Kannel and Bell (ref. 10) conducted a critical re- 
examination of the X-ray technique. Among the factors which may contribute to experi- 
mental inaccuracies, Kannel and Bell considered the influence of poor X-ray beam col- 
limation, the effects of possible surface reflections, and the consequences of lubricant 
X-ray absorption. They concluded that none of the aforesaid factors were of such sig- 
nificance as to seriously alter the accuracy of the X-ray measurements. 

Bell and Kannel (ref. 11) then examined several lubricant rheological factors which 
may help to account for the observed high-load effects. They considered in their analy- 
sis the effects of a non-Newtonian lubricant of the Ree-Eyring form whose viscosity ex- 
hibited a time -delay response to a sudden change of pressure. Their theory showed that 
for most practical cases, the time delay of pressure's effect on viscosity was suffi- 
ciently long as to not permit the lubricant's viscosity to change appreciably upon entering 
the contact zone. This resulted in a solution in which the film thickness exhibits a 
greater dependence upon load, h a than that predicted by all previous 

theories. Although Bell and Kannel' s solution is in better agreement with the observed 
trends, it still did not account for the rapid fall off of film thickness at the higher con- 
tact stress levels. 

Cheng (ref. 12) has attempted to extend present EHD theory to include the region of 
extremely heavy loads and high rolling speeds. Cheng determined that the introduction of 

ffPj+np-Pj) 

a composite exponential pressure-viscosity model (p = p Q e ) into the theory, 

as originally proposed by Allen, Townsend, and Zaretsky (ref. 13), in place of the con- 
ventional straight exponential pressure-viscosity model (p. = jj. Q e a ^) had negligible effect 
on the predicted film thickness. Furthermore, the introduction of a thermal reduction 
factor developed by Cheng in reference 14 to correct the isothermal theory for the in- 
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fluence of lubricant viscous heating at high-rolling speeds failed to significantly im- 
prove the agreement between theory and the experimental data within the heavy load 
regime . 

The variation of the lubricant's viscosity with changes in hydrodynamic pressure is 
known to play a vital role in the extent of surface separation within an EHD contact. Yet, 
there has not been enough work done toward establishing this relation under the appro- 
priate dynamic conditions. Even steady -state viscosity measurements have shown that 
not all lubricants obey the conventional straight exponential pressure -viscosity relation 
frequently used in classical EHD theory. A notable example of this can be found from 
the work of Chu and Cameron (ref. 15) in which it is shown that straight paraffinic 
mineral oils deviate substantially from the straight exponential relation (fig. 1). As an 
alternate, Chu and Cameron have developed a power law relation (p = p Q (l + Kp) n ) to fit 
the experimental pressure -viscosity data for the paraffinic based oils from reference 16. 

Dyson, Naylor, and Wilson (ref. 6) utilizing a similar power law pressure - 
viscosity relation developed a correction to be applied to the EHD theory of Grubin. 
However, a comparison between predicted film thicknesses with this correction and 
those obtained through capacitance measurement were not made due to uncertainties in 



Pressure above ambient, p, psi 


Figure 1. - Comparison of straight exponential pressure-viscosity relation with 
experimental data for paraffinic mineral oil (data from ref. 151. 


3 



estimating the value of the correction itself. Similarly, Cameron and Gohar (ref. 17) 
derived a generalized film thickness expression in which the power law pressure- 
viscosity relation could be used in place of the conventional straight exponential model. 

The objectives of the analysis reported herein were: (1) to incorporate into elasto- 
hydrodynamic theory a pressure -viscosity relation which obeys the power law model 
(fi = n Q (l + Kp) n ) developed in reference 15, and (2) to compare the numerical results 
from this analysis with the X-ray measurements of film thickness reported in refer- 
ence 7. 
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SYMBOLS 


minor semi-axis of Hertzian contact, m (in.) 
major semi-axis of Hertzian contact, m (in.) 
half width of Hertzian contact zone, m (in.) 

o 

modulus of 'elasticity of elements 1 and 2, N/m (psi) 



(psi) 


local film thickness, m (in.) 

film thickness in Cheng’s theory defined by eq. (29), m (in.) 

film thickness in Grubin's theory defined by eq. (18), m (in.) 

film thickness at inlet edge of Hertzian contact zone, m (in.) 

film thickness at point of maximum Hertzian pressure, m (in.) 

minimum film thickness, m (in.) 

film thickness with power law lubricant, m (in.) 

pressure-viscosity coefficient defined in eq. (3), m /N (psi ) 

constant defined indirectly in eq. (3) 

pressure -viscosity exponent defined in eq. (3) 

2 

local pressure above ambient, N/m (psi) 

2 

maximum Hertz stress, N/m (psi) 

2 

pressure at inlet edge of Hertzian contact zone, N/m (psi) 

critical pressure in composite exponential pressure -viscosity model, 
N/m^ (psi) 
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radius of elements 1 and 2 in rolling direction, m (in.) 



local temperature, K (°F) 
disk temperature, K (°F) 
l/2(Uj + Ug), m/sec (in. /sec) 

surface velocities of elements 1 and 2, m/sec (in. /sec) 

load parameter, w/E'R^ 

load per unit cylinder length, N/m (lbf/in.) 

dimensionless distance parameter, x/b^ 

distance from center of Hertzian contact zone, m (in.) 

2 2 

pressure -viscosity coefficient, m /N (in. /lbf) 

2 2 

equivalent pressure -viscosity coefficient defined in eq. (21), m /N (in. /lbf) 
function of T defined in eq. (4) 
constant defined indirectly in eq. (4) 

secondary pressure -viscosity coefficient in composite exponential pressure- 
viscosity model, m /N (in. /lbf) 

constant defined indirectly in eq. (4) 

-3 2 2 

local absolute viscosity, 10 N-sec/m or cP (lb-sec/in. or reyns) 

-3 2 2 

ambient absolute viscosity, 10 N-sec/m or cP (lb-sec/in. or reyns) 
Poisson's ratio of elements 1 and 2 
dimensionless pressure parameter, p^/Pjj z 


COMPARISON OF X-RAY DATA (REF. 7) WITH CONVENTIONAL 
ELASTOHYDRODYNAMIC THEORY 

Elastohydrodynamic film thickness measurements were obtained in reference 7 
using an X-ray rolling disk machine (fig. 2) with a synthetic paraffinic oil. The method 
of measuring film thickness with the X-ray technique comprises projecting X-rays be- 
tween the surfaces of the two contacting disks and detecting the amount of X-rays which 
are transmitted through the contact. Since the greatest constriction occurs at the trail- 
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^-Container-heater unit 

\ 



ing edge of the contact, the X-ray count thus becomes a measure of the lubricant's min- 
imum film thickness. 

The range of test conditions consist of disk temperatures from 339 to 589 K (150° to 
600° F), surface speeds from 9. 4 to 37.6 meters per second (370 to 1480 in. /sec) cor- 
responding to disk rolling speeds from 5000 to 20 000 rpm, and maximum Hertz stresses 
from 1.04x10^ to 2.42x10^ newtons per square meter (150 000 to 350 000 psi). Two 
crowned AISI M-50 steel disks, each with a radius of 1. 83 centimeters (0.72 in.) and a 
surface finish of 2.5x10"° to 5.0x10 centimeter (1 to 2 pin.) rms were used as the test 
specimens. Both crowned disks and crowned-cone disks (with a cone angle of 10°) were 
tested and no significant differences were found between the two sets of data. All data 
reported herein were generated with the latter disk configuration. The properties of the 
synthetic paraffinic oil which was used in the course of this study are listed in table I. 

A summary of the X-ray test results from reference 7 showing measured minimum 
film thickness h Q for a synthetic paraffinic oil plotted against maximum Hertz stress 
PHz for various mean surface speeds u and disk temperatures T q is presented in 
figure 3. The marked deviation of the X-ray data from theory at the higher stresses 
can be clearly seen from the dotted line in figure 3(a) whose slope represents the stress 
exponent of -0.22 from Cheng's theory (ref. 12). Evidently, within a heavily loaded con- 
tact, the film thickness dependence upon load as evidenced by the X-ray data is far 
greater than current EHD theory predicts. 

One possible explanation for this anomaly is that the lubricant's pressure varia- 
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Measured minimum film thickness, h-y pin. 








tion of viscosity does not follow the straight exponential model normally assumed by 
classical EHD theory. If the loss of pressure-viscosity dependence is to have an impact 
on the value of film thickness, it must occur at relatively low pressures. This is be- 
cause the nominal film thickness is basically established at the inlet of the contact zone 
where the hydrodynamic pressures within the lubricant are relatively low, generally, 
just a fraction of the maximum contact pressure. This deduction has been recognized 
by Bell and Kannel (ref. 11). It is also supported by Cheng's work (ref. 12) which 

showed that the loss of the lubricant's pressure -viscosity dependence at even moderate 

8 2 

pressures (i.e. , 4. 14x10 N/m (60 000 psi)) through the introduction of a composite ex- 

ap 1 +y(p-p 1 ) 

perimental pressure -viscosity model (p = p Q e ) did not produce a significant 

difference upon the calculated value of film thickness. 

The exact nature of the viscosity variation with pressure for the synthetic paraffinic 
oils is presently unknown. Lubricants are known to deviate from the familiar exponen- 
tial form: 


M 



( 1 ) 


where is the absolute viscosity at ambient pressure and p is the pressure above 
ambient. Chu and Cameron (ref. 15) reported that at high temperatures the straight 
paraffinic mineral oils do not obey the exponential relation shown in equation (1), even 
for low pressures as can be seen from figure 1. Based on the data from reference 16 
which contains pressure-viscosity information for a large number of fluids, Chu and 
Cameron developed a power law relation to fit the pressure viscosity data of the follow- 
ing form: 


M = m 0 (1 + Kp) n (2) 

where exponent n can be put equal to 16 and coefficient K is 

K = k(10)^ m 0 “°‘° 62 (3) 

and where 

0 = -% + XT) (4) 

2 -1 

For constant K in units of (N/m )~ , absolute viscosity in centipoises, and tem- 
perature T in K; 

k = 8.99X10" 9 


X = 0.0045 
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and constant /3 was taken to be -0. 75 for the straight paraffinic mineral oils. 

0 _i 

However, for constant K in units of psi , absolute viscosity ju Q in reyns, and 
temperature T in °F; 


k = 2.33X10" 5 


X = 0.0025 


and constant /3 Q becomes 0.4. 

It is of interest to examine the effects on film thickness of incorporating a power 
law viscosity variation with pressure of the kind described in equation (2) into an EHD 
inlet analysis of the Grubin type. 


POWER LAW PRESSURE-VISCOSITY RELATION 

Grubin (ref. 8) is credited with the first useful solution to the EHD problem of 
heavily loaded contacts. He had circumvented the lengthy numerical process attendant 
with the full EHD solution (e. g. , see ref. 9) through a clever simplification. Grubin 
assumed that the elastic displacements of the contacting lubricated surfaces under heavy 
loads will be essentially unchanged from the deformed shape in dry Hertzian contact. 
Grubin's geometrical model appears in figure 4. 
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With this assumption the film thickness distribution h(x) outside of the contact zone 
could be represented by the following Hertzian expression: 


h - h 


m 


- 2W 


r: 




X f X - 1 - In X 


+ Vx 2 - 1 j 


(5) 


where h m is the film thickness at dp/dx = 0, that is, center film thickness in this in- 
stance, and W, the nondimensional load parameter, and X are defined as follows: 


W = 


w 


and X = 


E'R' 


( 6 ) 


Here b^ is the semi-width of the contact zone and is equal to 


b h = 2R x wV2 


(7) 


and w is the load per unit cylinder length. Parameters R^ and E’ are defined as 


111 i 1 “ a i 1 ~ a 2 

— = — + — and — = + £. 


R i R 1 R 2 


E’ 7fE. 


nEc 


( 8 ) 


The hydrodynamic action of the lubricant within the contact zone is governed by the 
integrated form of the incompressible Reynolds equation as follows: 


^=12 M u( h “ hm 
dx 


(9) 


where dp/dx is the pressure gradient generated in the lubricant and u, the mean sur- 
face speed, is 


u 


u, + u„ 


2 


( 10 ) 


Equation (9) is developed on the assumption that the fluid is both Newtonian and incom- 
pressible. There is little to be lost by disregarding the effects of compressibility since 
the density increase normally experienced by the lubricant upon entering the contact 
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zone is relatively small and consequently has little influence on the ultimate separation 
of the contacting surfaces (see ref. 1) . 

Using Grubin's simplifying assumption, that the film thickness is constant over 
most of the contact zone, (i.e. , hj = h m , see fig. 4) and presuming that the lubricant’s 
viscosity is solely a function of pressure , that is 

p = p(p) (11) 

then equation (9) can be integrated from ambient pressure to the pressure at the inlet of 
the contact zone as follows: 



Utilizing the displacement expression shown in equation (5), Grubin numerically evalua- 
ted the integral appearing on the right-hand side of equation (12) for the practical range 
of h^ and found that 

W _ 3/2 x grubin's integral (13) 

x=-°° 



Dowson and Higginson (ref. 1) reported that Grubin's integral is fitted closely by 


0.0986 



- 11/8 


(14) 


therefore, equation (13) can be rewritten as 




0.1972 

R’ 


W-l/8/V\‘ 

w 


11/8 


(15) 
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and in terms of the inlet film thickness, rearranging equation (12) and combining with 
equation (15) yield 


h. 

1 


R' 


x 




- 8/11 


( 16 ) 


Equation (16) relates film thickness in terms of a generalized pressure -viscosity model. 
Grubin took the pressure variation of viscosity to be 

m(p) = P 0 e “ P (17) 


He also assumed that the pressure developed at the inlet was very large, that is, 
Pj^ — If these two conditions are applied, then equation (16) would reduce to the 
familiar Grubin formula 



1.87 


/ W' ' 

,-i/n 

«M 0 ( u i + u 2 >" 

\ E ' R x, 

) 

1 

r* 

« 

CQ 

1 


i 8/11 


(18) 


noting that h^ has been substituted for hu . Grubin's assumption that p i — 00 yields 
reasonable values of film thickness for lubricants whose viscosity increases rapidly 
with pressure, that is, those with large values of a. This is true analytically because 
the e -Qf P term which appears in the expression for film thickness when evaluating the 
pressure integral term in equation (16) is normally negligible for most practical situa- 
tions. (An alteration to Grubin's theory for a finite inlet pressure can be found in 
ref. 11.) However, for those lubricants which display a reduced viscosity -pressure 
dependence, the incorporation of a finite inlet pressure into the theory can have an 
appreciable effect on film thickness. Substituting the viscosity relation shown in equa- 
tion (2) into the pressure integral in equation (16) and performing the indicated integra- 
tion yield 


?1 dp _ 1 

fx(p) K(n - 1 V 0 

Substituting equation (19) into equation (16) results in a complete film thickness formula. 
Such a formula is 


(1 + Kp.) 


-n+1 


(19) 
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h / v 1/n 

i - 1 87 1 W ’ \ 

M 0 ( u i + u 2 } 

8/11 

> 

r 

i 

1- 1 


R' \ E'R' / 

x \ x/ 

2R i 

\ 

K(n - 1) 

(1 + Kp.) n-1 

J 


- 8 / 11 . 


(20) 


If the condition — 00 is applied to equation (20) and the result is then compared to 
Grubin’s formula in equation (18), a quantity, say a c , can be written as 

a c = K(n - 1) (21) 

where a c can be used in place of the customary a in Grubin's film thickness formula 
for lubricants obeying equation (2) in situations where the inlet pressure tends toward 
infinity. A similar expression was developed by Dyson, Naylor, and Wilson in refer- 
ence 6 and Cameron and Gohar in reference 17. 

For many practical EHD conditions, as will be shown later, the inlet pressure term 
in equation (20) can have an influence on the value of film thickness and should not be 
disregarded. For an approximate solution, it is convenient to express p^ as some 
fraction 4 of the maximum Hertz pressure. Thus, it is assumed that 

p i = ^Phz 

Now if exponent n in equation (2) is taken equal to its suggested value of 16 (ref. 15) 
and h is used to denote the film thickness for the power law lubricant, then equa- 

r 

tion (20) can be rewritten as 



before numerical comparisons can be made with the previous film thickness relation, 
estimates are needed for coefficients | and K. The value of £ is dependent upon the 
same parameters which influence film thickness, such as, surface speed, contact load, 
surface temperature, and type of lubricant as evidenced by the pressure distributions 
presented in reference 18. Nevertheless, taking 4 = 0.3, as an approximation, should 
yield informative results representative of the conditions under consideration. A simi- 
lar approach was taken in reference 11. 
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For purposes of an approximate solution, a reasonable estimate for coefficient K 
can be made as follows: At low temperatures, less than about 100° F (311 K), refer- 
ence 15 (fig. 1) shows that the viscosity predicted by the power law model in equation (2) 

does not deviate appreciably from the value predicted by the straight exponential model 

8 2 

in equation (1) for pressures to 80 000 pounds per square inch (5.5x10 N/m ). Hence, 
at these conditions, 


(1 + Kp) 16 * e ap (24) 

or 

(1 + Kp) « e"^ 16 (25) 

Expanding the right side of equation (25) in a Taylor series, we obtain 

e op/16 = 1 + op + ( ap/16 ) 2 + + (ofp/16) n ( 26 ) 

16 2'. ’ nl 

Comparison of equations (25) and (26) shows that, for low pressures, a good first-order 
approximation for K is 


K « — (27) 

16 

Currently, the only available published data for a of a synthetic paraffinic oil 

Q 

appears in reference 19 where a was reported to be 1.3x10 square meter per newton 
(0.92X10 -4 psi -1 ) at 298 K (77° F). Thus solving equations (3), (4), and (27) simul- 
taneously at T = 298 K (77° F) one finds that 

/3 = -0.48 for SI units (28) 


or 


/3 = 0.67 for U.S. customary units 

Using this value of /3 Q together with equations (3) and (4) , one can calculate the value 
of K at the various disk temperatures as given in table II. 
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RESULTS AND DISCUSSION 


A comparison of the power law film thickness formula developed in equation (23) 
with the experimental X-ray data (ref. 7) is presented in figures 5 to 7 and tabulated in 
table in. To contrast the present results with previous EHD theory, the isothermal film 
thickness equation of Cheng (ref. 20) has been included in these plots. Cheng's film 
thickness equation for line contact or more precisely elliptical contact with b/a > 5 
can be written 



°' 74/ Ph z V 0 - 22 

E’ 


(29) 


Cheng's theory (ref. 20) 

Power law theory 

X-r ay data (from fig. 3) 
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Figure 5. - Comparison of predicted minimum film thickness 
with X-ray data for synthetic paraffinic oil at disk temperature 
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T 0 = 339 K (15(P F). 
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In computing Cheng's minimum film thickness and the minimum film thickness from the 
present analysis, a correction factor of 0.8 has been applied to adjust the predicted 
center film thickness to minimum film thickness. The geometry of a typical EHD con- 
tact showing the relative size of the inlet, central, and minimum film thicknesses, is 
presented in figure 8. 



Figure 8. - Geometry of typical elastohydrodynamic contact 
showing relative size of inlet (h j), central (h c ), and mini- 
mum (h 0 ) film thicknesses. 


Before one can make a fair comparison of Cheng's formula with the other data, an 
estimate of the pressure -viscosity coefficient a at various temperatures is required. 
Unfortunately, data of this kind is not currently available for synthetic paraffinic oils. 
To overcome this obstacle, a first-order approximation of the variation of a for the 
synthetic paraffinic oil was modeled after the trend exhibited by the paraffinic mineral 
oils determined by the ASME Research Committee on Lubrication (ref. 16). The 
pressure -viscosity coefficients at ambient pressure for five paraffinic mineral oils 
labeled 31G to 35G were determined at temperatures of 273, 298, 311, 372, and 491 K 
(32°, 77°, 100°, 210°, and 425° F) . Figure 9 shows a plot of a at the aforementioned 
temperatures normalized by a at 298 K (77° F) as a function of temperatures for the 
five oils. By applying this factor (a /a 2 g 8 K or oi/oirjrj o F ) to the synthetic paraffinic 
oil's pressure -viscosity coefficient a curve of a against temperatures as presented in 
figure 10 can be obtained. A comparison of the exponential pres sure -viscosity model 
(eq. (1)) with the power law pressure-viscosity model (eq. (2)) utilizing the estimated 
pressure -viscosity coefficients a and K appears in figure 11. 

Incorporating the estimated value of a interpolated from figure 10 into Cheng's 
equation leads to a reasonably good trend and magnitude agreement with the X-ray data 
at moderate contact pressures, p Hz = 1.04x10® newtons per square meter (150 000 psi), 
at least for the higher temperatures as can be seen from figures 6 and 7. This agree- 


17 



Pressure-viscosity coefficient, o, psi 






CSJ 




Pressure above ambient, p, psi 


Figure 11. - Comparison of power law and straight exponential vis- 
cosity variation with pressure for synthetic paraffinic oil. 


ment would not have been possible had the reduction of a with increasing temperature 
not been considered. Although the slopes of Cheng’s film thickness formula are com- 
patible with the X-ray results toward the lower load regime, they deviate substantially 
as the load increases. 

The present theory from equation (23) differs little from previous theories at the 
lower temperatures. As previously discussed, the lubricant's pressure variation of 
viscosity can be adequately represented by either the exponential or power law models 
to moderate pressures at low disk temperatures as evidenced by figure'll. At the 
higher disk temperatures, the present theory displays a slightly stronger load depend- 
ence than previous isothermal EHD theories, but it also fails to account for the rapid 
reduction of film thickness at the higher loads. 

The physical mechanism responsible for the high load dependence displayed by the 
measured minimum film thickness remains unexplained. There is some evidence (e.g. , 
ref. 19 for the case of point contact) which suggests that the minimum film thickness 
ma y be inherently more sensitive to applied load than is the central film thickness. If 
this is the case, then a simple inlet-type EHD analysis would not be capable of reflecting 
this increased film thickness load sensitivity. Nevertheless, it is evident that the use 
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of current EHD film thickness formulas, for most design purposes, will yield unreal- 
istically high values of the minimum film thickness at maximum Hertz stress much 

Q 

above 1.04x10 newtons per square meter (150 000 psi). 


SUMMARY OF RESULTS 

An isothermal elastohydrodynamic (EHD) inlet analysis of the Grubin type which 
considers a power law viscosity variation with pressure (p = /i q (1 + Kp) n ) and a finite 
pressure at the inlet edge of the Hertzian contact zone was performed. Results from 
the present EHD inlet analysis were compared to the X-ray test data reported in NASA 
TN D-6411 and to the results from a conventional EHD analysis for line contact. Both 
present and previous EHD theories were numerically evaluated with estimated values 
of pressure-viscosity coefficients a and K for a synthetic paraffinic oil. The results 
of this comparison are as follows: 

1. Analytical results from the elastohydrodynamic theory utilizing a power law 
pressure -viscosity relation display a slightly stronger film thickness dependence upon 
load than previous isothermal theories which assume a straight exponential pressure- 
viscosity relation. However, the model did not adequately define the decrease in film 
thickness with contact stress displayed by the X-ray test data. 

2. The minimum film thickness predicted by conventional elastohydrodynamic theory 

exhibits reasonably good trend and magnitude agreement with the test data at moderate 

9 2 

contact stresses (1.04x10 N/m (150 000 psi)) for the higher disk temperatures (above 
422 K (300° F)) . 

3. For contact stresses above 1.04x10^ newtons per square meter (150 000 psi), the 
conventional elastohydrodynamic theory fails to account for the radical reduction of film 
thickness with increasing load displayed by the test data. 

Lewis Research Center, 

National Aeronautics and Space Administration, 
and 

U.S. Army Air Mobility R&D Laboratory, 

Cleveland, Ohio, September 19, 1972, 

501-24. 
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TABLE I. - PROPERTIES OF THE SYNTHETIC PARAFFINIC OIL 


-fi 9 

Kinematic viscosity, cS (or 10" m /sec), at - 


233 K (-40° F) >100 000 

311 K (100° F) 443 

372 K (210° F) 40 

478 K (400° F) 5.8 

589 K (600° F) a 2.3 

Flash point, K (°F) 541 (515) 

Fire point, K (°F) 589 (600) 

Autoignition temperature, K (°F) 703 (805) 

Volatility (6.5 hr at 533 K (500° F)), wt. % 14.2 

Specific heat at 533 K (500° F), J/(kg)(K) (Btu/(lb)(°F)) 2910 (0.695) 

Thermal conductivity at 533 K (500° F), J/(m)(sec)(K) (Btu/(hr)(ft)(°F)) 0.12 (0.070) 

Specific gravity at 533 K (500° F) 0.71 

Pressure-viscosity coefficient at 298 K (77° F) , m 2 /N (psi“*) ^1 . 3x10"® (0. 92x10’^) 


a Extrapolated. 
^From ref. 19. 


TABLE II. - ESTIMATED VALUES OF 
PRESSURE-VISCOSITY COEFFICIENT 
K FOR A SYNTHETIC 
PARAFFINIC OIL 


Disk 

temperature 

Pressure-viscosity 

coefficient, 

K 

K 

°F 

m 2 /N 

psi" 1 

339 

150 

60.7X10" 11 

41.9X10" 7 

366 

200 

48.5 

33.4 

422 

300 

29.6 

20.4 

478 

400 

17.5 

12.1 

505 

450 

13.3 

9.2 
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TABLE III. - COMPARISON OF PREDICTED MINIMUM FILM THICKNESS with X-RAY DATA 
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Crowned-cone disk data from ref. 































































TABLE III. - Continued. COMPARISON OF PREDICTED MINIMUM FILM THICKNESS WITH X-RAY DATA 
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Crowned -cone disk data from ref. 




TABLE in. - Continued. COMPARISON OF PREDICTED MINIMUM FILM THICKNESS WITH X-RAY DATA 
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Crowned-cone disk data from ref. 7 






































































TABLE III. - Concluded. COMPARISON OF PREDICTED MINIMUM FILM THICKNESS WITH X-RAY DATA 



NASA-Langley, 1973 


is E-6973 
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Crowned -cone disk data from ref. 
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